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Abstract: Leading-twist distribution amplitudes (DAs) of light mesons like tt, p etc. de- 
scribe the leading nonperturbative hadronic contributions to exclusive QCD reactions at 
large energy transfer, for instance electromagnetic form factors. They also enter B decay 
amplitudes described in QCD factorisation, in particular nonleptonic two-body decays. 
Being nonperturbative quantities, DAs cannot be calculated from first principles, but have 
to be described by models. Most models for DAs rely on a fixed order conformal ex- 
pansion, which is strictly valid for large factorisation scales, but not always sufficient in 
phenomenological applications. We derive models for DAs that are valid to all orders in the 
conformal expansion and characterised by a small number of parameters which are related 
to experimental observables. 
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1. Introduction 

Light-cone distribution amplitudes (DAs) appear in the description of hard exclusive QCD 
processes by factorisation, schematically 



where the labels i and j refer to the hadrons in the incoming and outgoing states, re- 
spectively, (j)(n) is the DA that describes the probability amplitude for a hadron to be 
found in the parton state n and T(Q 2 ;ni,rij) is the perturbative function that describes 
the hard scattering between the partons; Q 2 is the hard momentum transfer. The symbol 
(g) indicates a convolution, i.e. a sum or integral over the parton degrees of freedom that 
correspond to the states n« and rij. Factorisation implies that the process can be split into 
two separate regimes, the hard-scattering subprocess T, which involves a short-distance 
momentum transfer and can be calculated in perturbation theory, and the long-distance 
hadronisation of the partons emerging from the short-distance process, described by the 
process-independent DAs (p. The expansion of the amplitude in ( |L1| ) is ordered by inverse 
powers of the momentum transfer Q 2 , which corresponds to a light-cone expansion of the 
process in terms of contributions of increasing twist. In this paper we are concerned with 
the leading- twist DA of light mesons. Although we shall focus on the pseudoscalar mesons 
7T and K, the models we suggest are equally well applicable to — and in fact have already 
Q been applied to — vector mesons (p, uj, K* and (p). 

The leading-twist meson DA of a meson is related to its Bethe-Salpeter wave function 
by integrating out the dependence on the transversal momentum k±, 



amplitude ~ JJ 4> oxl t,j{nj) <8> T(Q 2 ; 



(1.1) 
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where u is the longitudinal momentum fraction carried by the quark (and u = 1 — u that 
carried by the antiquark) . Originally introduced in the context of hadron electromagnetic 
form factors or the pseudoscalar-photon transition form factor B], light-meson DAs have, 
in recent years, attracted increasing interest also in B physics due to their appearance in 
QCD sum rules on the light-cone |l], [6| and factorisation formulas for B decay ampli- 
tudes Q and form factors ||. The theory of meson DAs is well understood |9|, |ll| and 
suggests their parameterisation in terms of a partial wave expansion in conformal spin. One 
advantage of this expansion is that the contributions of higher conformal spins to the con- 
volution integral and hence the physical amplitude are suppressed by the highly oscillating 
behaviour of the corresponding partial waves, another one that conformal symmetry en- 
sures, to leading-logarithmic accuracy, multiplicative renormalisation of each partial wave. 
The combination of both features suggests the construction of models for DAs based on 
a truncated conformal expansion, where only the first few waves are included, typically 
one to three. Not much is known, in general, about the amplitudes of these partial waves. 
For the ir and 77 some experimental information on the tt (17)77* transition form factor is 
available Jl^], which one can use to extract values of the first two amplitudes. From the 
theoretical side, there exist a few dated lattice calculations for the second moment of the 



7r DA p3l , and a recent retry [14]; unfortunately, these results are still preliminary and 
cannot yet be used in phenomenological applications. Other theoretical calculations, for 
both pseudoscalar ||, 10, |l5|, 16] and vector mesons [jH], [lTf], use the method of QCD sum 
rules, which turns out to be not very suitable for higher moments; we will come back to that 
point in Sec. I). 1 It is probably fair to say that at present some information is available on 
the lowest moments of the n DA, but much less so for other mesons. In view of the growing 
demands on B physics to deliver "precision results", cf. e.g. Ref. fl20fl , and the prominent 
role of QCD factorisation in testing the mechanism(s) of flavour violation, it then appears 
timely to assess the actual theoretical uncertainty of hadronic decay amplitudes induced by 
the truncation of the conformal expansion and to devise alternative models for light-meson 
DAs that do not rely on conformal expansion, but establish a closer connection between 
the characteristics of DAs accessible in "classical" applications of pQCD in hard exclusive 
processes and their use in B physics. This is precisely what we aim to achieve in this paper. 

Our models for the leading-twist DAs of the tt and K are based on the fall-off behaviour 
of the nth Gegenbauer moment of these amplititudes, a n , in n, which we assume to be 
power-like. We shall argue that such a behaviour can be justified for instance from the 
known perturbative contributions to the a n . We formulate our models in terms of a few 
parameters, notably the first inverse moment of the DA, which, for the tt, is directly related 
to experimental data, and the strength of the fall-off of the a n . The models can be summed 
to all orders in the conformal expansion and give predictions for the full DA at a certain 
scale. Depending on the parameters, the models also predict a nonstandard behaviour of 
the DAs at the endpoints, which affects for instance the scaling, in m&, of the B— > light 
meson form factors at zero momentum transfer. 

Our paper does not aim to give a fully-fledged analysis incorporating all available 

1 These problems have been evaded in Refs. Jl8|, |l^] in a modifed version of QCD sum rules using nonlocal 
condensates. 
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constraints on the it DA from low-energy experimental data as scrutinised in Refs. [19, 21 



22]. Rather, we aim to work out the gross features of leading-twist DAs, which are likely 
to apply to all light pseudoscalar and vector mesons. 

Our paper is organised as follows: in Sec. [2] we define the leading-twist DA of pseu- 
doscalar mesons and discuss its partial-wave expansion in conformal spin. In Sec. |3| we 
derive models for the DA based on the fall-off behaviour of higher partial-waves and for- 
mulate constraints on the model parameters. In Sec. |] we investigate the dependence of one 
important quantity in B physics, the semileptonic form factor /+~ > ' r (0), on the model DAs, 
and in Sec. || we study their effect on CP asymmetries and branching ratios in nonleptonic 
B decays. We summarise and conclude in Sec. ||. 

2. Definitions and Conformal Expansion 

Distribution amplitudes are defined in terms of matrix elements of nonlocal operators near 
the light-cone. For pseudoscalar mesons P in particular one has 



b P (u) + - m 2 pX 2 A P (u) + 0(x 4 ) 



(0\q 1 (x)^ 5 [x,-x]q 2 (-x)\P) = if P p^ / due l ^ x 

Jo 

I^K-xJ due*** [M P (u) + 0(x 2 )] , (2.1) 
I px J 

where £ = 2u — 1 and we use [x, y] to denote the Wilson- line connecting the points x and 



[x,y] = Pexp 



igf dt(x-y)^(tx + (l-t)y) 
Jo 



(2.2) 



In Eq. ( |2.1| ), (f)p is the leading twist- 2 DA, whereas Ap and Bp contain contributions from 
higher-twist operators. The corresponding definitions of vector meson DAs can be found 



in Ref. g0|. 

The extraction of the leading behaviour of the matrix elements on the light-cone yields 
ultraviolet divergences, whose regularisation generates a nontrivial scale-dependence that 
can be described by renormalisation group methods || [|. Conformal invariance of QCD 
allows one to express the DA in terms of a partial wave expansion, also called conformal 
expansion, in contributions from different conformal spin, which do not mix with each other 
under a change of scale. This is true to leading logarithmic accuracy, but no longer the case 
at higher order, as the underlying symmetry is anomalous. For the leading-twist DA 4>(u), 
for both pseudoscalars and vectors, the conformal expansion is in terms of Gegenbauer 
polynomials C^!" 1 ', 

oo 

0(n 5/ u 2 ) = 6u(l-u)^2a n ( f i 2 )C^ 2 (2u-l). (2.3) 

n=0 

The coefficients a n , the so-called Gegenbauer moments, renormalise multiplicatively to 
leading logarithmic accuracy: 

MQ2) = W) (2 . 4) 
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with 0o = 11 — (2/3)nj. The one-loop anomalous dimension is given by [23] 



7 w = 8CF ^ ( „ +2)+7E _^_ (n + i) ^ +2) ). (2 . 5) 

For 7r, p, cj, i], rf and 0, G-parity ensures that a Q( jd = and that the DA is symmetric 
under u <-> 1 — u, whereas for if and if* the nonzero values of a Q dd induce an antisymmetric 
component of the DA. ao = 1 is fixed by normalisation, 

p(u,fi 2 ) = 1, 

whereas all other a n are intrinsically nonperturbative quantities. As they do not mix under 
renormalisation, Eq. ( t2.3|) is well suited to construct models for eft: truncating the series 
after the first few terms, typically three, yields a parameterisation of the DA that is "stable" 
under a change of scale, except for the numerical values of a n . Despite there being no small 
expansion parameter in the game, such a truncated conformal expansion is often - but not 
always - a meaningful approximation to the full DA, as we shall see below. 

As the anomalous dimensions are positive (except for 7q ' 1 = 0), the contributions of 
higher a n get suppressed for large scales, and for fi 2 — > oo the DA approaches the so-called 
asymptotic DA 

M («) = 6u(l - u). (2.6) 

For many processes involving DAs, in particular in B physics, it is usually argued that 
a truncated conformal expansion be sufficient for the calculation of physical amplitudes as 
long as the perturbative scattering amplitude is "smooth" — the reason being the highly 
oscillatory behaviour of higher order Gegenbauer polynomials. It is actually instructive to 
quantify this statement, for example for the simplest, but phenomenologically very relevant 
case of one meson in the initial or final state, so that the convolution integral reads 

1=1 du<p{u)T(u) } (2.7) 
J o 

where T is the perturbative scattering amplitude. 
We distinguish the following cases: 

(i) T is nonsingular for u € [0,1]; 

(ii) T has an integrable singularity at one of the endpoints; 

(iii) T contains a nonintegrable singularity at one of the endpoints. 
As a typical example for case (i), consider T(u) = y/u, which yields 

f 1 , , ( s (-l)"36(n + l)(n + 2) 

Jo dU(p{UNu ~ (2n - l)(2n + l)(2n + 3)(2n + 5)(2n + 7) ° n ' 

This result implies a strong fall-off ~ 1/n 3 of the coefficients of higher Gegenbauer moments 
a n : assuming dj = 1 for all i, already the first three terms in the sum account for 98.8% of 
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the full amplitude. In reality the convergence will be even better as all existing evidence 
points at \a n \ <C 1 for n > 1. 

As an example for case (ii) choose T(u) = In u, which yields 



The singularity at u = evidently worsens the convergence of the series; again assuming 
a% = 1, the first three terms now overshoot the true result by 35%. In order to approximate 
the full amplitude to within 5% one now has to include nine terms, but the convergence 
will again be better in practice, thanks to the fall-off of a n in n. 

Case (hi) is more complicated and depends on the asymptotic behaviour of the a n . For 
T(u) = 1/u, for instance, one obtains 



Here the amplitude is finite only if the a n fall off sufficiently fast in n. For stronger endpoint 
divergences the coefficients multiplying a n start to grow in n and for T ~ 1 /u 2 the integral 
diverges, even for the asymptotic DA, which would indicate a breakdown of factorisation 
for that process. 

The conclusion to be drawn from this discussion is that models for eft, based on a 
conformal expansion that is truncated after the first few terms, are indeed appropriate for 
cases (i) and (ii), but less so for case (in). Convolutions with T ~ 1/u are actually very 
relevant both in hard perturbative QCD reactions, e.g. 77* — > tt, and in B physics, e.g. B — > 
7T7T (7j. Given the different weight the a n do have in different convoluted amplitudes, and 
the fact that their impact is highest in convolutions of type (iii), it appears not unreasonable 
to base a parameterisation of <fi not on, say, a2 and 04, setting all a n >4 = , but rather on 
the "worst case scenario" of Eq. fl2.8| ), where all a n enter with the highest possible weight 
factor. This is the basic idea of our models of leading-twist DAs we shall elaborate on in 
the next section. 

3. Models for <j>„ and (fi K 

There are basically two properties that any viable model for the leading-twist tt DA <p n 
must fulfill: 

(a) in the limit of large energies, fj? — > 00, the DA must approach the asymptotic DA 
4>ir(u, [J 2 = 00) = 6uu Q; 

(b) the first inverse moment, du^(u)/u, which is related to the 7177* transition form 
factor, must exist. 

Condition (a) must evidently be fulfilled also for all other light-meson DAs; condition (b) 
should be fulfilled in general if QCD factorisation in B decays @ is to make sense. Obviously 
any model that is based on a truncated conformal expansions fulfills both constraints, and 





(2.8) 
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in addition predicts that for u — > 0, 1, <p ~ u(l — independent of the factorisation scale. 
We shall show that this prediction is, in general, not fulfilled for models that are not 
truncated at fixed order in the conformal expansion. 

The starting point for our models is the relation between the first inverse moment of 
4> and the sum of all Gegenbauer moments: 

As mentioned before, one has to distinguish between mesons for which a dd vanish due 
to G-parity, such as the ir, and strange mesons, for which the odd moments induce an 
antisymmetric part of the DA. It is the former case we shall study first. Available exper- 



imental data for the tt, as summarised in Ref. 22], point at a value of A around 1.1 at 
the scale ji ~ 1.2 GeV, which implies that the infinite sum be convergent. Hence, even 
at the comparatively low scale 1.2 GeV the am must fall off fast enough in n. Assuming 
an asymptotic equal-sign behaviour for large n, a power- like fall-off as a2n ~ l/n p with p 
slightly larger than 1 is one of the slowest possible fall-offs. It turns out that DAs defined 
by a power-like fall-off of the Gegenbauer moments can actually be summed explicitly: 
using the generating function of Gegenbauer polynomials, 

_. oo 

/K - t)= (l-2 P + t ^ = g Cy(g '" 



the DA with moments 



for n even, 



is given by 



n (n/b + l) a, 

= jf dt (- ln *) a_1 {f( 2u - 1» tl/b ) + f( 2u - !> - tl/b : 



In the same way one obtains a DA with alternating-sign behaviour of the Gegenbauer 
moments, 

(_l)«/2 



(n/b + l) a ' 



for n even, 



as 



= fSj£ dt (~ ln {f( 2u - itVb ) + f( 2u ~ !' ~ itllb ) 

The corresponding values of A are A+ fo = (b/2) a ((a, 6/2) and A~ 6 = (6/4) a (C(a, 6/4) - 
C(a, 1/2 + 6/4)), where C( a i s ) = Ylh=o + s ) a ^ s ^ ne Hurwitz ^ function. In order to 
obtain models for arbitrary values of A, we split off the asymptotic DA and write 

0± b (A) = 6uu + A ± ~ 1 - 6uu) , valid for a > 1 and 6 > 0. (3.2) 

a,b 

Is there any reason why the a n should fall off as inverse powers in n — other than that 
the corresponding DAs can be expressed in closed form via Eq. ( [3.2P ? The answer to this 
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question relies on the behaviour of a n (/j, 2 ) under a change of scale. Evidently the models 
are defined at one particular scale, for instance the hadronic scale [i ~ 1.2 GeV. At larger 



scales Q , the a n change according to Eq. (2.4). For large n, we have 



(n) 

7o 



8C F lnn + 0(l) 



and 



a n (Q 2 



1 



■a n (fj, 2 



n iC F //3 ln(l/L) 

with L = a s (Q 2 )/a s (n 2 ). That is: perturbative leading-order scaling induces a power-like 
fall-off of the a n , at least for large n. We take this as indication that such a behaviour is 
indeed intrinsic to QCD. One more consequence of scaling is that for Q 2 —* oo, i.e. L — ► 0, 
the suppression of higher a n is power-like with a power that approaches infinity, so that 



(Q 2 — > oo) = (f> (a —* oo) = 6u(l — u). 



Hence 4> ab as defined in (3.2) approaches the asymptotic DA in this limit, which implies 
that condition (a) follows from (b). 

Eq. (3.2) implies that the asymptotic DA is recovered for A = 1, and also from (j)^ b 
in the limit a — ► 1. The models are valid only for a > 1, as otherwise A^ fc diverges, 
or, equivalently, (f> ab does not vanish at the endpoints u = 0, 1. In Fig. |l] we plot several 



examples of (f> a b for a fixed value of A; it is evident that the two models 



a,b 



and • 



a,b 



have a 



rather dissimilar functional dependence on u; in particular 4> h turns out to be nonanalytic 



wc 



at u = 1/2 for a < 3. The dependence of 4> a b on b is illustrated in Fig. |2|. In Fig. 
show the possible values of the lowest Gegenbauer moments a2,4 that can be obtained for 
fixed values of A and 6, but different a. For 



one always has 02 < A — 1, for 



(Z2 > A — 1; in the limit a —> 00 both branches meet and one obtains the truncated NLO 
conformal expansion with 02 = A — 1 and 04 = 0. 
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Figure 1: Left: Examples for model DAs <p^ b as functions of u, for a = 1.5,2,3,4, b — 2 and 



1.2 (solid curves), as compared to the asymptotic DA (dashed curve). For a 



approaches the asymptotic DA. Right: the same for (f> a 2 . Note that for a < 3, <p a b is nonanalytic 
at u = 1/2 and displays a pronounced "spike". 



Before discussing constraints on the model parameters A, a and b, let us shortly 
comment on the behaviour of the DAs near the endpoints u = 0, 1. In many applications of 
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Figure 2: Left: Examples for model DAs (j)^ b as functions of u, for b — 0.2,2,20 and A = 1.2 
(solid curves), as compared to the asymptotic DA (dashed curve). For b — + oo, <j)^ b approaches the 
asymptotic DA, as A^^ — > oo. Right: the same for 4>^ b . 
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Figure 3: Values of and 0,4 for different 
parameters A and a of (f> a 2 . The curves in 
the upper right and lower left quadrant cor- 
respond to 4i 2 , those in the lower right and 
upper left one to </>~ 2 . 
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Figure 4: Comparison of <pf 2 i u ) (solid 
curves) with the asymptotic DA (dashed 
curve) for A = 1.04. The "spike" at u = 1/2 
characteristic for <p~ causes these models to 
significantly deviate from the asymptotic DA 
even for A close to 1. 



DAs it is assumed - implicitly or explicitly - that (f) ~ u(l— u) near the endpoints. Although 
this is evidently the case for any model based on a truncated conformal expansion, it does 
not apply to our models, at least not for all values of a. A closer inspection shows that 
4> + ~ yj u(l — u) for a = 2 and ~ u(l — u) ln(n(l — u)) for a = 3. It is only for a > 3 that 
the DAs behave in the "canonical" linear way. Are there are any rigorous arguments why 
the DAs should behave linearly near the endpoints, even for low hadronic scales? As far as 
we could trace the origin of the argument in favour of linear behaviour, it was stated first 
in Ref. U and relies on the fact that QCD sum rules for the moments of the DA exhibit 
the following behaviour: 



em 
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Two comments are in order here. First, an exact calculation for arbitary </> shows that 



so that the correct statement is that (£ n ) ~ 1/n 2 implies that 4>'(u) exists at the endpoints 
- which in turn indeed implies a linear behaviour, since 0(0) = = 0(1) as long as A < oo. 
Second, the conclusion of the authors of Q is based on results obtained from the leading- 
order perturbative contributions to QCD sum rules. NLO expressions have been obtained 
in |l6|, [TtJ and yield ~ u(l — u) ln 2 (u/(l — u)), which upsets the linear approach to the 
endpoints and is equivalent to a n ~ 1/n 3 . The large- n behaviour of nonperturbative terms 
cannot be obtained from these sum rules, 2 but we see no reason why it should not follow 
the a n ~ 1/n 3 behaviour of NLO perturbation theory or even introduce 1/n 2 scaling. We 
hence conclude that the standard assumption of ~ u(l — u) is not rigorously justified at 
hadronic scales and that all corresponding conclusions, in particular the scaling behaviour 
of the B — > light meson form factors at zero momentum transfer, /(0) ~ l/ m b > first 
derived in p4| ], should be taken cum grano salis. We also would like to add one more 
remark about the expectation that for small A — 1 the DA should be "very close" to the 
asymptotic DA. In Fig. || we show that this statement is indeed true for + models, but 
not for 0~, which are characterised by their nonanalytic behaviour at u = 1/2. Of course 
the asymptotic DA is recovered for A — 1 —* 0, but not in an analytic way. 

Let us now discuss the experimen- 
tal restrictions on A, a and b for the 
7r. First, we fix the reference scale 
as u = 1.2 GeV. As the data are too 
scarce to constrain all parameters in 
a meaningful way, we fix b = 2. As for 
A, we require it to be larger or equal 
1. This is the case if a<i > 0, which is 

indeed the common overall conclusion of all determinations available in the literature, as 



A+ 



6 



oo 



2.04 
1.04 



1.58 
1.11 



1.43 
1.16 



1.36 
1.19 



1.33 
1.22 



1.27 
1.27 



Table 1: Upper bound on A , for various values of 
a, as implied by 0^(1/2) > 0.9. The reference scale is 
H w 1.2 GeV. 



compiled in Ref. [22]. An upper bound on A can be inferred from experimental data or 
from requiring a,2 < 0.2, which is again the maximum value allowed by most analyses. 3 
For + the restriction on 02 implies A < 1.2, which coincides with the allowed range of A 
found in |]l9|| , whereas for the resulting range of A is smaller. One more constraint on 
(j)^ is the allowed range of </>„-( 1/2), which follows from light-cone sum rules for the ttNN 
coupling fll5f : 

0.9 < &r (1/2,1 GeV) < 1.5. 

For A > 1, 0(1/2) is always smaller than 1.5, so it is only the lower bound that is relevant. 
In Tab. |l| we list the corresponding maximum A for various a. Evidently the constraints on 

2 The reason being that a n for large n are intrinsically nonlocal quantities, which cannot be obtained 
from a local operator product expansion (OPE). The ineligibility of the local OPE manifests itself as 
contributions to a n that scale as positive powers in n, which is incompatible with a finite value of A. 

3 With the notable exception of Chernyak and Zhitnitsky || , whose results are however by now generally 
considered to be excluded by experiment. 



- 9 - 



A + are weaker than the ones discussed before, but for <p~ the minimum value of 4> n (l/2) 
poses a nontrivial constraint on A. 

It is possible to further refine the constraints for the ir as it was done in e.g. [19, 21]. 
In this paper, however, we are not so much interested in the ir, but rather in the gross 
characteristics of the DAs, which are likely to be valid also for other pseudoscalar and 
vector mesons. We hence refrain from pursuing that line of investigation, but summarise 
the main constraints on the symmetric parts of our model DAs, which are likely to be valid 
also for other mesons: 

• 1 < A < 1.2 for < a>2 < 0.2 for <p^ 2 : this is based on the observation Q that DAs 
of mesons with higher mass tend to become narrower; 

. 1 < Min(1.2, A max ) for 0" 2 , with A max given in Tab. § 

• 6 = 2, lacking further data. 

Let us now turn to 4>k, the twist-2 DA of the K meson. It differs from </>„- by the fact 
that now also odd Gegenbauer moments contribute. Models for the antisymmetric part of 
the DA can be constructed in a similar way as before as 4 



The models are characterised by c and yield a\ = (2/3) c . Models for the antisymmetric 
part of 4> with arbitrary a\ can then be defined as 

^ : = o 1 (3/2)V±(«). (3.4) 

Examples for such models are shown in Fig. || 

Numerical values for have been discussed in ||; we quote af (1.2GeV) ~ 0.15. For 
(j>K, the total Atot can then be written as 

^tot,± _ ^ ^_ ^asym,± 

where A = ^ a e ven is the contribution of the symmetric part of the DA to A tot and A asym,± 
is defined as 

A asym ' + = f du'^P- = - ai (3/2) c C(c,3/2), 
Jo 3u 

A asym <~ = ^ du = - ai (3/4) c {C(c, 3/4) - C(c, 5/4)} . 

Jo 3ti 



fciu) = ^ J^dti-lnty- 1 (/(2n-l,^)-/(2u-l,- v / t)), 

CM = f^y T^-lntr 1 (/(2«-l,iV^)-/(2«-l,-iVt)). (3.3) 



4 In complete analogy to the symmetric part of the DA we could introduce one more parameter d that 
would correspond to b. In view of the near complete absence of any information on antisymmetric DAs, we 
refrain from writing down the corresponding formulas and set d = 2 from the very beginning. 
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Figure 5: Models for the antisymmetric contributions to the twist-2 DA for a\ = 0.15. Left: ip+ 
as function of u for c € {2, 3, 5} (solid curves), dashed curve: ip^. Right: tp~ as function of u for 
c e {1,2,3,5} (solid curves), dashed curve: ip^. Like the symmetric models <f>~ , ip~ is nonanalytic 
at u = 1/2 for c < 3. 

4. Results for f* ^(0) 

One important application of factorisation in B physics is the calculation of the weak decay 
form factor ff^ 71 from QCD sum rules on the light-cone. It is beyond the scope of this 
paper to review the method of QCD sum rules on the light-cone (LCSRs), for which we 
refer to Ref. pi]. Instead, we would like to stress that LCSRs offer a means to calculate 
nonperturbative hadronic quantities like form factors within a controlled approximation, 
which relies on the factorisation of an unphysical correlation function, whose imaginary 
part is related to the hadronic quantity in question. 
The B — * ir form factors are defined as 



2 2 

(p + Pb)^ B 2 n q„ 



/rV) + ^y^/o^V), (4-1) 



{<K{p)\V,\B(jp B )) 

with q = pb — p- Within the LCSR method f+~o* can be related to a correlation function 

n^ps) = i J d A x^ x {K{p)\TV^x)j B (m (4.2) 

= U + (q 2 ,p B )(p + p B ) fl + U^(q 2 ,p 2 B )q l , , 

where jb = mbdvy$b is the interpolating field for the B meson. For unphysical p 2 B <C m 2 , 
H± can be expanded around the light-cone as 

Ul c (q 2 ,p 2 B ) = V f 1 duT^\u,q 2 ,p 2 B ,^ n \u^), (4.3) 

where the sum runs over contributions of increasing twist and the term in n = 2 is just 
the leading-twist contribution. J 1 ^ 2 ' 3 ^ are known to 0(a s ) ||, whereas T± is known at 
tree- level |5[. We have extended the calculation performed in || to include contributions 
up to aio- The LCSR for the form factor / + depends on the spectral density /?+ c of 11^. 
in p 2 B and reads 

r s o 

e ~m 2 B /M 2 <■ "/;— - > / 



m 2 B f B f^(q 2 ) = I" dse- s l M2 p\ C {s,q% (4.4) 

J ml 
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where M 2 , the so-called Borel parameter, and sq, the continuum-threshold, are sum rule 
specific parameters. Using the central values of these parameters as obtained in ||, M 2 = 
9.2 GeV 2 and s = 33.9 GeV 2 , and the model twist-2 DAs 0± 2 we obtain the values of 
jB— >tt(q) s ] lown [ n Yig. |6|. The dependence of >7r (0) on A, with a fixed, is linear, as the 




A 



Figure 6: /f^O) calculated fr om LCSRs, using the twist-2 DAs (f>^ 2 (solid lines) and <j) a2 
(dashed lines). Left: /+ _ " r (0) as function of A for a — 3 fixed; right: ff^^(0) as function of a for 
A = 1.1 fixed. The endpoints of the dashed curves are set by the constraint 0(1/2) > 0.9. 

values of the individual a n >2 are just rescaled by a common factor. If, on the other hand, 
A is kept fixed and a is being varied, higher a n are increasingly suppressed with increasing 
a, so that for a — > oo the form factor obtained using cp + approaches that obtained using 
, as is clearly visible in Fig. |(| 
Figure || shows that the theoretical uncertainty induced by A and a is about ±5%. 
This has to be compared with the final value of f+ >7r (0) stated in Ref. ||: 0.258 ± 0.031, 
i.e. a 12% uncertainty, which includes also the variation of other input parameters of the 
LCSR. 



5. B — > nn and B — > Ktt in QCD Factorisation 

Another important application of leading-twist DAs are nonleptonic B decays treated in 
QCD factorisation [0]. Recent experimental data point at a failure of QCD factorisation to 
explain the observed branching ratios and CP asymmetries in B — > tttt and B — > irK, which 



has prompted a number of authors to explain this effect by new physics, e.g. Ref. [25|, but 
has also motivated other authors to investigate the impact of nonfactorisable corrections 
to QCD factorisation, which are suppressed by inverse powers of nib p6fl - In this section 
we investigate the effect of nonstandard DAs on the predictions of QCD factorisation. The 
factorisation formulas depend on twist-2 DAs of the ir and K, but also on the form factor 
/+"^ 7r (0); f° r which we use the results from LCSRs, cf. Sec. 

Let us first study the time-dependent CP-asymmetry in B — > 7r + 7r~, which is defined 

as 

Agp = r (fo->7r+7r ) - r(fl -> 7T+7T ) = 
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where our interest is in the mixing-induced asymmetry S nn which depends on the unitarity 
triangle (UT) angles (3 and 7 via 

_ 2ImA^ 2l0 e~^ + P/T 

a ™-l + \ X ^\2> A ™~ e dh + p/T' [ ' 

In the limit where the penguin-to-tree ratio P/T is zero this reduces to S n7T = sin 2a. 
Although P/T is highly suppressed, it is not negligible and can be expressed, in QCD fac- 
torisation, in terms of the CKM parameters and the Gegenbauer moments a n . Neglecting 
the small contributions from weak annihilation terms 0, we find that P/T is given by a 
ratio of polynomials in a n with complex coefficients. For the asymptotic DA, for instance, 
one has 

P _ 1 (-1.36 - 0.318) - (1.46 + 0.37i) a + (0.22 + 0.0H) 

T ~~R~ b (13.0 - 0.08i) + (15.5 + 0.03i) a - (0.56 - 0.15i) a? Q ' ^ ' ' 



where R b = V ' P 2 + V 2 = 0.40 @. 

The variation of S W7T in terms of a and A is shown in Fig. [/]. For 1 < A < 1.2 as 
suggested in Sec. ^, the variation is of order 1% and becomes more significant only for 
unrealistically large values of A; the convergence for A = 1 corresponds to the asymptotic 
wave function. The current experimental results [28] are 

= -0.30 ±0.17 ±0.03 (BaBar), 5™ = -1.00 ± 0.21 ± 0.07 (Belle), (5.4) 

which could only be accommodated using very extreme values of A. For example, taking 
a model with same-sign fall-off, to reproduce the BaBar result with a = 2 would require 
a value of A > 10 to be within the la band and A = 20 to approach the central value. 
Even for higher values of a, a minimum A ~ 16 is needed to approach S nn = —0.3, and 
would produce Gegenbauer moments that are significantly outside the known constraints, 
for example with a = 6 and A = 16, we find a2(2.2 GeV) = 7.7 and 04(2.2 GeV) = 0.5. 
Obtaining a value of SVtt = —1-0 is not possible for values of A > 1, which is required to 
keep 02 positive. 

The effect of the DA model on the branching ratios of B — > tt + it~~ is significantly more 
pronounced than for the CP-asymmetry. The central value of the branching ratio for the 
asymptotic DA is 

BR(B -» tt+O = 5.5 x 10- 6 |0.25e i l5 ° + e~^\ 2 , (5.5) 

where 7 = 60° ± 7° , and the explicit dependence of the branching ratio on the Gegen- 
bauer moments is again a polynomial in a n : 

BR(B vr+vO » 1.50 x 10~ 6 |1 + (1.17 - 0.01*) a - (0.06 + O.Oli) a§ (5.6) 
+ (1.16 - O.Oli) a 2 - (0.02 + O.Oli) a\ - (0.08 + 0.02t) a a 2 + ... | 2 , 

where the terms of first order in a n come from nonspectator interaction and those of second 
order from the hard-gluon exchange. Figure shows the variation of the branching ratio 
for the two model DAs, which has to be compared to the current experimental world 
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average BR(B — > 7r + 7r~) = (4.6 ± 0.4) x [29]. We see that large increases from the 

asymptotic value are possible by increasing A, with values within the physical range, with 
a considerable effect of around 30% possible with the alternating-sign fall-off model. On 
the other hand, the effect of a nonasymptotic DA is to increase the branching ratio, and 
hence leads the result of QCD factorisation even further away from the experimental result. 



A c p (B->JT + JT ) 
-0.56 




BR (B->7T + 7T ) (l(T b ) 
9 




Figure 7: Left: Acp(tttt) as a function of A for a = 2, 3 ... 6, shown for models DAs with same-sign 



fall-off (dashed lines) and alternating sign fall-off (solid lines). Right: the same for BR(_B 
Both sets of curves converge for increasing values of a. 



7T r 7T 
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The situation for B — * irK decays is complicated by the presence of the K DA (f>K, 
for which very little is known about the Gegenbauer moments. As a result of this, we use 
the same moments for as in (f) w with the addition of the parameter af , for which we 
use the value of af(1.2 GeV) = 0.15, as discussed in Sec. ||. Concentrating on the decay 
Bq —* ir + K~ and its CP-conjugate, we consider first the direct CP asymmetry, recently 
reported in Ref. J30| as 



A 



ttK 
CP 

attK 
A CP 



-0.133 ± 0.030 ± 0.009 (BaBar), 
-0.101 ±0.025 ±0.005 (Belle). 



(5.7) 



In the framework of QCD factorisation, Acp can be written in terms of real, strong in- 
teraction parameters (derived from the factorisation coefficients) and pure CKM variables. 
We can neglect the annihilation contributions, but the electroweak penguins play a crucial 
role and must therefore be included, which yields 



attK 

a cp 



tan 2 9 C Rb (sin 7) r(irK) 
1 + tan 2 9 C R b (cos 7) ? (ttK) ' 



(5. 



where 9 C is the Cabibbo angle and r, r'(irK) contain the QCD effects which depend on the 
DA model parameters A and a. 

Figure |8| shows the dependence of the direct CP-asymmetry in Bq — > tt + K~ on A 
and the fall-off parameter a. The first point to note is that Aqp differs in sign from 
the experimental value (|5.7| ), which is in agreement with the findings of [0]. The precise 
value of Aqp does depend on A, but it is impossible, even for extreme DAs, to obtain the 
experimentally observed negative sign. This is also emphasised by plotting the dependence 
of the asymmetry on the UT angle 7, as expressed in (|5.8| ), which is also shown Fig. || 
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A CP (B->7T + fC) 





30 ° 60 ° 90 ° 120 ° 150 ° 180 ° 



Figure 8: Left: the variation of as a function of A for a = 2, 3 ... 6, for same-sign fall-off. 
Right: the dependence of on the UT angle 7 for curves at a = 2, A — 1 (lowest curve) to 
A = 4. 



There is a 10% increase in the value of Acp(ttK) between the asymptotic form with A = 1 
and A = 2, more significant changes only occurring outside the physical range of A. 

As in the B — > it + it~ case, the effect of our model DAs on the branching ratio of 
-Bo ~~ > it K~ is much more pronounced than for the CP-asymmetry and can exhibit up 
to 30% change from the asymptotic value within the physical ranges of a and A. The 
branching ratio is given as 



BR(B ttK) oc (F B " 



)) 2 im;i 2 



1 2 



(5.9) 



with 6km = tan# c 2 i?ft and c,c'{ttK) are derived from factorisation coefficients, containing 
all the dependence on hadronic parameters and the DAs <f> n k- Using the asymptotic DA 
we find BR(i3 — > irK) = 13.58 x 10 -6 , and the variation from this when higher-order 
Gegenbauer moments are included is shown in Fig. ^. The experimental world average, 



BR(B->7T*fT) (10 ) 
20r 



BR (B->7r + K~ ) (10" 

21r 
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Figure 9: BR(B — * n + K~) as a function of A for a — 2, 3 ... 6, for models DAs with same-sign 
fall-off (left plot) and alternating sign fall-off (right plot). The experimental average is also marked 
for comparison. 

also shown in Fig. g, is (18.2 ± 0.8) x 10~ 6 There can be significant changes to the 

asymptotic value, especially in the model with alternating sign fall-off. The experimental 
average is shown for comparison and can be accommodated with A ~ 1.8 for the same sign 
fall-off, or A ~ 1.4 with alternating sign fall-off. 
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In Fig. [T(J we also show the depen- 
dence of the direct CP-asymmetry in B —* 
ttK on ai, the Gegenbauer moment that 
parameterises the antisymmetric part of 
the K DA, whose actual value is around 

0. 15. It is evident that one can shift Acp 
into the "right" direction by decreasing 

01, but again extreme values of a\ would 
be needed to obtain a negative CP-asymmetry. 

Our conclusion from this investiga- 
tion is that the discrepancy between the 
experimental values of BRs and CP-asym- 
metries and their predicted values in QCD 
factorisation can not be attributed to the 

uncertainties in the leading-twist DAs: for B — * 7r + 7r~ it is impossible to reproduce the 
data, whereas for the direct CP-asymmetry in B — > irK highly unrealistic values of A and 
a would be needed in order to reconcile the theory predictions with experimental data. 




Figure 10: The direct CP-asymmetry in Bq 
tt + K~ as function of a\ for a = 3 and A 
{1,2,3,4} (from bottom to top). 7 = 60°. 



6. Summary and Conclusions 

We have presented models for the symmetric part of leading-twist light-cone distribution 
amplitudes (DAs) of light mesons which depend on three parameters. Two of these pa- 
rameters control the fall-off behaviour of the Gegenbauer moments a n in n, whereas the 
third one, A, is given by the first inverse moment of the DA and parameterises the max- 
imum possible impact of higher Gegenbauer moments on the actual physical amplitude 
described in factorisation. We have also developed similar models for the antisymmetric 
part of the DA, which is relevant for K and K* mesons; these models are normalised to 
d, the first Gegenbauer moment. For the ir DA, for which experimental data exist, we 
have formulated constraints on the model parameters which are likely to be valid also for 
other meson DAs. We have argued that these models are better suited to estimate the true 
hadronic uncertainty of processes calculated in factorisation than the standard truncated 
conformal expansion and have studied these uncertainties for two quantities, the B — > tt 
weak decay form factor '"'(O) and the CP asymmetry in B — > -kit and B — > Kit. For the 
former, the theoretical uncertainty induced by the model-dependence of the DA is smaller 
than that due to other parameters and approximations. For nonleptonic decays calculated 
in QCD factorisation we find that the branching ratios are more sensitive to the precise 
values of the model parameters than the CP-asymmetries. In both decays channels it is 
however impossible to explain the experimental data by nonstandard DAs, which indicates 
the presence of nonnegligible nonfactorisable contributions — a conclusion that agrees with 



the findings of other authors, e.g. Refs. |2g, |33| . 



Our models should prove particularly useful for describing DAs of mesons other than 
the 7T which are also symmetric by virtue of G-parity, but for which no experimental or 
other reliable theoretical information is available — in particular the p, uj and (j). For 
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these particles, we argue that existing theoretical indications from local QCD sum rules 
g, 0, point at the DAs being narrower than that of the it, which implies the allowed 
values of A being smaller than 1.2 at the scale 1.2 GeV. The same results also imply 02 to 
be positive, which entails A > 1. For the parameter a, which controls the fall-off of the 
a n in n, we have found that the perturbative contributions to QCD sum rules indicate it 
to be 3 [11], but that also smaller values of a are not excluded unless one can rigorously 
prove that the leading-twist DAs must behave as ~ u(l — u) near the endpoints u = 0, 1, 
at all scales, which would imply a > 4. One more relevant constraint for the models <fi~ b 
with an alternating-sign fall-off of the Gegenbauer moments comes from the requirement 
07r(l/2) > 0.9 at the scale fj, = 1 GeV. A positive value of this quantity is also required by 
QCD sum rules on the light-cone for the couplings gDD*-n [0] and qdd*p Iffil , and hence 
also very likely to be the case for other mesons. But even without these constraints being 
taken literally, our models provide a way to test the impact of nonasymptotic DAs on 
physical amplitudes without recourse to conformal expansion. 
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